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. ABSTRACT 

The masses of the central supermassive black holes and the Eddington ratios for a sam- 

■ pie of 65 young radio galaxies [27 gigahertz-peaked spectrum (GPS) and 38 compact steep- 
I spectrum (CSS) sources] are estimated by various methods. We find that the average BH mass 

of these young radio galaxies is <log i\/bh>— 8.3, which is less than that of radio loud QSOs 
and low redshift radio galaxies (<logAfbh>— 9.0). The CSS/GPS sources have relatively 
high Eddington ratios with average ratio <log Lboi/iEdd>= —0.56, which are similar to 
those of narrow line Seyfert 1 galaxies (NLSls). It suggests that the CSS/GPS sources may 
(-H ' not only be in the early stage of radio activities, but also in the early stage of their accretion 

" ' activities. 

We find that the young radio galaxies as a class deviate systematically from the il/bh — cr* 
O ' relation defined by nearby inactive galaxies, when using (J[q mj as a surrogate for the stellar 

^ , velocity dispersion cr* . There is no significant correlation between the deviation of the [O III] 

^ ■ emission line width. Act = a^Q mj — criprcd]' ^nd the jet/accretion power, where cr [pj-cd] 

calculated from the Tremaine et al. relation using the estimated BH masses. However, we find 
that the deviation Act in young radio galaxies is well correlated with the Eddington ratio, and 
this correlation is found to be similar to that of radio quiet AGN (QSOs/NLSls) where the 
\ radio jet is absent or weak. We suggest that the accretion activities may still play an important 

■ role in shaping the kinematics of [O III] narrow line in these young radio galaxies. 

' Key words: accretion, accretion discs - ISM: jets and outflows - galaxies: active - radio 

I continuum: galaxies - quasars: emission lines 
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1 INTRODUCTION frustration scenario, the radio source is confined and frustrated by 

a cocoon of dense gas and dust for its entire life. However, the far- 
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Gieahertz-peaked spectrum (GPS) and compact steep-spectrum j n j ■ r j ^ 

" . ^ . ^ n^ - inirared fluxes and near-inirared propertities ot CSS sources are 

(CSS) radio sources constitute a large fraction (~ 407o) of the pow- 



, a-sur TT -1 ' f roughly consistent with those of radio sources with extended jets, 

5-H erful(loBPi4 J5 10 W Hz ) radio source population. Their ra- ... , , , ■ ■, i ■ i ■ 

rrt ' which suggest that they have similar dust content in the circumnu- 

. . I dio spectra are simple and convex with peaks close to 1 GHz and , . ^ , _ . . nr,nr, j -ir ■ . i ir>no^ 

clear environment (e.g., Fanti et al. 2000; de Vries et al. 1998), 



100 MHz for GPS and CSS sources respectively. Morphologically, 
the GPS sources are generally less than 1 kpc in projected linear ^ the evolutionary hypothesis of the youth scenario is cor- 
size and CSS sources are less than 20 kpc (see O'Dea 1998 for ^ct, smdy of the early-phase of the radio activity and feedback in 
a review). Two main theories have been proposed to explain the CSS/GPS sources may impact our understanding of galaxy evo- 
CSS/GPS phenomenon: the youth scenario (Fanti et al. 1995) and lutio"- There is much evidence to suggest that, in a significant 
ths frustration scenario (van Breugel 1984). The age of CSS/GPS fraction of the radio source population, the radio sources are trig- 
sources has been estimated using both the lobe proper motion (e.g., gered by merger events of two or more galaxies (e.g., Heckman 
Owsianik & Conway 1998) and the radiative ages (Murgia et al. et al. 1986). There are many compact radio sources exhibiting 
1999). The ages derived are less than a few thousand years, there- features attributed to mergers, such as double nuclei, tidal tails, 
fore strongly supporting the youth scenario in the interpretation arcs of emission and distorted isophotes, suggesting that these 
of their compactness. In fact, it has been suggested that compact sources are possibly observed relatively shortly after the merger 
CSS/GPS sources perhaps eventually evolve into extended, edge- event (Stanghellini et al. 1993). The merger of galaxies will deposit 
brightened Fanaroff-Riley 11 sources (Snellen et al. 2003). In the large quantities of gas and dust into the nuclear region, the cen- 
tral black hole (BH) will grow rapidly through the merger-induced 
accretion, and eventually form strong winds/outflows during fast 
* E-mail: qwwu@shao.ac.cn accretion (e.g., Di Matteo et al. 2005). Narrow line Seyfert 1 galax- 
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ies (NLSl) are a particular class of active galactic nuclei (AGN). 
The most widely accepted paradigm for NLSls is that they accrete 
at close to the Eddington rate and have smaller BH masses com- 
pared broad line AGN, which might be in the early stage of AGN 
evolution (e.g., Mineshige et al. 2000; Boroson 2002; Grupe 2004; 
Bian & Zhao 2004). Most of the radio loud NLSls galaxies are 
compact, steep-spectrum sources, which are similar to CSS/GPS 
sources (e.g., Gallo et al. 2006; Komossa et al. 2006). Therefore, 
compact radio sources (CSS/GPS) may not only be in the early 
stage of radio activity, but also in the early stage of BH accretion. 

The mass of supermassive BH in the AGN is a key param- 
eter to understand the nuclear energy generation as well as BH 
formation and evolution. There are several methods for estimating 
the BH masses in AGN. In quiescent galaxies, dynamical model- 
ing of either stellar kinematics (e.g., Kormendy 1988) or gas mo- 
tions (e.g.. Harms 1994) is used to determine central BH masses 
(e.g., Kormendy & Gebhardt 2001 for a review). The correlation 
between the bulge stellar velocity dispersion (j» and the central 
BH mass Mbh is now well established in nearby inactive galax- 
ies (Mbh — o", relation, e.g., Gebhardt et al. 2000a; Fenarese & 
Merritt 2000; Tremaine et al. 2002). Interestingly, the Afbh - cr, 
relation for normal galaxies also extends to AGN (e.g., Gebhardt 
et al. 2000b). Unfortunately, the bulge stellar velocity dispersion in 
bright AGN is generally difficult to measure directly. In order to 
derive the BH mass in a larger sample of AGN, the width of the 
[O III] line emitted by the narrow-line region is frequently used 
as a surrogate for the stellar velocity dispersion, since the kine- 
matics of [O III] emission line is believed to be dominated by 
the gravitational potential of the bulge. Nelson & Whittle (1996) 
found o-[o III] — cr. in a sample of 66 Seyfert galaxies, where 
f^[o III] = FWHM(0 III)/2.35 and FWHM is full width at half- 
maximum of the emission line. Nelson (2000) and Boroson (2003) 
found that Afbh and g^q ni] are strongly correlated and consistent 
with the Tremaine et al. (2002) relation, although with substan- 
tial scatter. The reverberation mapping (RM) method provides a 
virial estimate of the BH mass in type 1 AGN. Using this tech- 
nique, the size of the broad line region (BLR) can be measured 
using the time lag between the variations of the continuum and 
emission line fluxes. The BH mass can then be estimated with 
Mbh ~ iiBLRVVG, if we assume the eas near the BH is viri- 
alized, where i?BLR is the BLR size, V is the characteristic veloc- 
ity (e.g., FWHM), and G is the gravitational constant (e.g.. Bland- 
ford & McKee 1982; Peterson 1993; Peterson 2004, and references 
therein). The most important result of the RM method is the dis- 
covery of a simple scaling relationship between the BLR size and 
the optical luminosity, which provides a useful secondary method 
of mass determination (e.g., Wandel et al. 1999; Kaspi et al. 2000; 
Vestergaard 2002; Kaspi et al. 2005; Greene & Ho 2005a; Kong 
et al. 2006). This technique has been widely used in the literature 
(e.g., Gu et al. 2001; Wu et al. 2004). 

Despite the widespread use of cr[o m] as a surrogate of a* in 
AGN, there are also several exceptions. In particular, one should 
be cautious where the [O III] kinematics may be affected by other 
drivers (e.g., outflows/winds) apart from the first driver of grav- 
itational potential of bulge. For example, it is found that NLSls 
as a class deviate systematically from the Mbh — f* relation seen 
in normal galaxies, when assuming cr* = ayo m] (e.g., Grupe & 
Mathur 2004; Bian & Zhao 2004). However, Barth et al. (2005) 
measured ct, in 19 NLSls, and found that these NLSls still follow 
the Mbh — relation. The width of low-ionization lines (e.g., [S 
II]A6716, 6731) or the core component of [O III] (after removal 
of asynmietric blue wings) of NLSls is still a good proxy for stel- 



lar velocity dispersion (e.g., Greene & Ho 2005b; Komossa et al. 
2008). Therefore, the accretion activity (e.g., accretion-driven out- 
flows/winds) may affect the kinematics of the [O III] emission line 
(e.g., Greene & Ho 2005b; Bian et al. 2005; Komossa et al. 2008). 
The other possible exception is the [O III] narrow emission line in 
young radio galaxies (CSS/GPS sources) and some luminous linear 
radio sources, which, normally, have highly broadened line profiles 
(e.g., Gelderman & Whittle 1994; Nelson & Whittle 1996; Tad- 
hunter et al. 2001; O'Dea et al. 2002; HoU et al. 2003, 2006, 2008). 
High resolution optical and radio imaging show that the [O III] op- 
tical narrow emission line is strongly aligned with each other in 
CSS sources at all redshifts (e.g., de Vries et al. 1997; de Vries et 
al. 1999; Axon et al. 2000; Privon et al. 2008). The expanding ra- 
dio jets provide a convenient driving mechanism for the observed 
nuclear outflows which may broaden the [O III] emission line (jet- 
driven outflows, e.g., Tadhunter et al. 2001; Holt et al. 2003, 2006, 
2008). Labiano (2008) also found that the [O III] luminosity is well 
correlated with the CSS/GPS source size, which seems to support 
the jet-driven scenario. The ionization diagnostics on different nar- 
row lines found that mix of shock and photo-ionization give the 
best explainations in the extended emission line gas (e.g., Inskip 
et al. 2002; Labiano et al. 2005; Holt et al. 2006; Humphrey et al. 
2008). 

Much of the work on the compact radio galaxies has concen- 
trated on the radio wavebands and optical emission lines, but rela- 
tively little attention has been paid to their central engine properties 
(e.g., BH masses and accretion modes). In particular, the accretion- 
driven outflows/winds may also affect the kinematics of the [O III] 
emission line in the radio loud CSS/GPS sources in addition to 
jet-driven outflows, if their Eddington ratios are high as those of 
radio quiet NLSls. Therefore, it is important to test these two sce- 
narios in young radio galaxies. In this paper, we present the BH 
masses, Eddington ratios, and the possible relations between the 
[O III] kinematics and the accretion/jet properties for a sample of 
CSS/GPS sources. Throughout this paper, we assume the following 
cosmology: Ho = 70 km s^^Mpc^^, fio ~ 0.3 and = 0-7. 



2 SAMPLE 

Our goal is to investigate the nuclear properties (e.g., accretion 
and jet) and their possible relation to the kinematics of the [O III] 
emission line in the compact radio galaxies (CSS/GPS sources), 
whose selection is unbiased with respect to the nuclear properties 
themselves. For each source, we must be able to estimate the BH 
mass, and/or bolometric luminosity, jet power, the [O III] kinemat- 
ics (e.g., FWHM). We present 38 CSS sources, of which 27 sources 
are selected from Spencer et al. (1989) and the other 1 1 sources are 
selected from different papers (see Table I for references). We also 
selected 27 GPS sources from the literature from various surveys 
(e.g., Snellen et al. 1999, 2003), of which two sources (OQ 172 and 
4C 12.50) were defined as CSS sources in Spencer et al. (1989), but 
later were defined to be GPS sources according to their radio spec- 
tra and sizes (e.g., O'Dea 1998 and references therein). Therefore, 
our sample includes 65 compact radio sources (38 CSS -I- 27 GPS 
objects). We note that our sample is a characteristic, rather than 
complete sample of compact radio sources. The histogram distri- 
butions of the redshift for our sample are shown in Figure I. The 
dotted-line is for CSS sources, the dashed-line is for GPS sources, 
and the solid-line is their sum. 

Table I lists our sample with the relevant information. 
Columns (l)-(3) give the object's name, redshift, and references 
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Table 1. BH masses, Eddington ratios, [O III] kinematics, and jet powers. 



Source 


z" 
(2) 


Refs. 
(3) 


A4h 

lo2 (Mr^') 

(4) 


J-bol 
-^Edd 

lo? 
(5) 


Method l*" 
(6) 


Method T 
(7) 


Refs.'* 
(8) 


""[O HI] 
log (km/s) 
(9) 


Refs. 
(10) 


Qjct 

los? (orff /s) 

(11) 


Band 

ri2) 












CSS sample 














3C 43 


1.459 


S89 


9.2 


-1.0 


Lopt+V 


L5100 


M06 






47.20 


160M 


3C 48 


0.369 


S89 


8.8 


-0.30 


Lnp+V{B.f3) 


L-Ril 


I91,C97 


2.84 


G94 


46.49 


178M 


3C 67 


0.310 


S89 


8.1 


-1.52 


Lopt+V 


L5100 


S07 


2.41 


G94 


45.71 


160M 


3C 93.1 


0.244 


S89 


7.1 






Fitting 


W02 


2.42 


G94 


45.49 


160M 


3C 138 


0.759 


S89 


8.7 


-0.25 


Lopt+V 


Fitting 


W02 


2.72 


G94 


46.71 


160M 


3C 147 


0.545 


S89 


8.7 


-0.90 


Lh/3+V^(H/3) 


L-Rp 


L96,M06 


2.71 


G94 


46.85 


178M 


3C 186 


1.063 


S89 


8.9 


-0.30 


Lh;3+V(H/3) 


Lnp 


H03 


2.67 


H03 


46.88 


151M 


3C 190 


1.195 


S89 


7.8 


-0.13 


Lnp+V{B.P) 


Lnp 


H03 


2.61 


H03 


47.07 


190M 


3C 191 


1.956 


S89 


9.7 


-1.27 


iciv+V(CIV) 


Lciv 


M06,C91 






47.49 


160M 


3C 213.1 


0.194 


S89 


9.1 


-0.42 


A4'bh--f'ion 


T ■ *^ 


D03 


2.52 


G94 


45.07 


160M 


3C216 


0.670 


S89 


7.0 


-0.30 


Lh;3+V(H/3) 


L-Rp 


L96 






46.63 


178M 


3C241 


1.617 


S89 


7.8 


0.56 


Lh/3+V(H/3) 


Lup 


H03 


2.63 


HOB 


47.71 


151M 


3C 268.3 


0.372* 


S89 


7.8 








NED 


2.46 


G94 


45.87 


178M 


3C 277.1 


0.320 


S89 


7.6 


-0.14 


Lh/3+V^(H/3) 


Ln/s 


G94 


2.34 


G94 


45.65 


178M 


3C 286 


0.849 


S89 


8.5 


-0.09 


Lopt+V 


L5100 


M06 


2.31 


G94 


46.87 


151M 


3C 287 


1.055 


S89 


9.6 


-0.85 


Lopt+V 


L5100 


M06 






46.87 


151M 


3C 293 


0.045 


S04 


8.0 


-1.78 




L5100 


W02,X08 






44.22 


151M 


3C 298 


1.439 


S89 


10.4 


-0.83 


Lopt+V 


L5100 


M04 






47.68 


178M 


3C 303.1 


0.270* 


S89 


8.4 








NED 


2.54 


G94 


45.54 


151M 


3C 309.1 


0.905 


S89 


9.1 


-0.45 


Ly_p+V(B.I3) 


Lhp 


L96 






46.95 


178M 


3C 343 


0.988 


S89 


7.5 


0.14 


Lopt+V 


Fitting 


W02 






46.81 


178M 


3C 346 


0.162 


S89 


8.8 


-2.54 




L5100 


NED,M04 


2.60 


G94 


45.16 


160M 


3C 380 


0.692 


S89 


9.4 


-1.15 


Lopt+V 


L5100 


M04 


2.41 


G94 


47.07 


178M 


3C 454 


1.757 


S89 


8.6 


-0.58 


LciY+V{ClW) 


Lciv 


B90 






47.25 


160M 


3C 455 


0.543 


S89 


8.5 


-1.62 


Lopt+V 


L5100 


M04 


2.25 


G94 


45.84 


160M 


3C 459 


0.220 


H08 


8.5 




Mr-M^^ 




NED 


2.63 


H08 


45.79 


160M 


2342+821 


0.735 


S89 


7.5 


0.18 


Lopt+V 


Fitting 


W02 


2.39 


G94 


46.00 


151M 


4C 14.82 


0.235 


S89 


7.9 


-0.79 


LH0+y(H/3) 


Lhi3 


C97 


2.33 


G94 


45.20 


160M 


PKS 0252-71 


0.563* 


H08 


8.0 






NED 


2.48 


H08 


46.11 


408M 


PKS 1151-348 


0.258 


M97 


9.1 


-1.54 


Lopt+V 


Fitting 


W02 






45.53 


160M 


PKS 1221-423 


0.171 


J05 


7.8 




Mr-M^^ 




NED 


2.53 


J05 


45.04 


160M 




1 .Do / 




Q 1 


n Qi 


Mfl-iUbh 


T s 1 no 








Al 1 1 


1 DUiVl 


PKS 1549-79 


0.152 


H08 


8.0 


-0.32 


Mfl-Mbh 


LhI3 


H06 


2.77 


H06 


45.29 


408M 


PKS 2004-447 


0.240 


G06 


7.6 


-0.26 


T, t+V 

^opt^y 


Fitting 


W02 






44.94 


408M 


SDSS 1722+565 


0.425 


K06 


7.5 


-0.14 


Lopt+V 


L5100 


K06 


2.32 


K06 


44.48 


330M 


RXJ0134-4258 


0.237 


K06 


7.6 


0.01 


Lopt+V 


L5100 


K06 


2.34 


K06 


44.73 


4850M 


TeX 11111+329 


0.189 


K06 


7.3 


-0.21 


Lopt+V 


L5100 


K06 


2.80 


K06 






IRAS 2018-2244 


0.185 


K06 


6.5 


1.25 


Lopt+V 


L5100 


K06 


2.47 


K06 


44.0 


352M 












GPS sample 














B 0019-0001 


0.305 


S03 


8.4 




Mbh-f^* 




S03 




... 


45.09 


178M 


B 0830+5813 


0.094 


S03 


7.6 




Mbh-fT* 




S03 






42.84 


325M 


B 0941-0805 


0.228 


S03 


7.7 




Mbh-T* 




S03 










B 1819+6707 


0.221 


S03 


8.4 




Mbh-o-* 




S03 


<2.2 


S99 


44.25 


325M 


B 1144+3517 


0.063 


S03 


8.2 




Mbh-fT* 




S03 






42.93 


151M 


B 1946+7048 


0.101 


S03 


7.9 




Mbh-o-* 




S03 


<2.15 


S99 


43.56 


325M 


B 2352+4933 


0.238 


S03 


8.2 




Mbh-o-* 




S03 






44.65 


151M 


B 0537+6444 


2.417 


S99 


1.1 


-0.29 


iciv+V(CIV) 


LciY 


S99 






45.39 


325M 


B 0544+5847 


2.860 


S99 


8.0 


-0.91 


iciv+V(CIV) 


Lciv 


S99 






45.77 


325M 


B 0601+5753 


1.840 


S99 


7.9 


-0.53 


Lciw+V{Cl\) 


Lciv 


S99 






45.16 


325M 


B 0755+6354 


3.005 


S99 


8.3 


-0.41 


LciY+V{ClY) 


Lciv 


S99 


... 




45.72 


325M 


B 0758+5929 


1.977 


S99 


8.1 


-1.05 


iciv+V(CIV) 


Lciv 


S99 






45.84 


325M 


I) nQTA_i_'7n/i ^ 
£5 Uozo+/U4j 


Z.UUo 


coo 


c ^ 

O.J 




i^ciY+y (S'^y ) 


Lciv 








4 J .40 




B 1538+5920 


3.880 


S99 


7.7 


-0.16 


iciv+V(ClV) 


Lciv 


S99 






45.99 


325M 


B 1550+5815 


1.320 


S99 


8.9 


-0.65 


LMgn+V(MgII) 


i'Mgll 


S99 






45.42 


325M 


B 1622+6630 


0.202 


S99 


7.7 


-0.94 


Lhc«+V(H«) 


Lfla 


S99 


<2.15 


S99 


43.26 


325M 


B 1642+6701 


1.895 


S99 


9.1 


-0.79 


iciv+V(CIV) 


Lciv 


S99 






45.89 


325M 


B 1647+6225 


2.180 


S99 


7.2 


-0.38 


iciv+V(CIV) 


Lciv 


S99 






45.50 


325M 


B 1746+6921 


1.886 


S99 


7.5 


-0.22 


Lciv+V(CIV) 


Lciv 


S99 






45.65 


325M 
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Table 1. (Continued) 



Source 


z 


Rets. 


AJbh 


-^Edd 

log 


Method 1 


Metnoa z 


Rets. 


o"[0 III] 


Rets. 


Vjot 


Band 








log (Mq) 








log (km/s) 




log (erg/s) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


B 1945+6024 


2.700 


S99 


7.5 


-0.70 


iciv+V(CIV) 


^'CIV 


S99 






45.61 


325M 


B 1958+6158 


1.820 


S99 


7.8 


-0.21 


iciv+V(CIV) 


Lew 


S99 






45.53 


325M 


B 0237-233 


2.223 


S99 


9.5 


-0.77 


Lciv+V{CIY) 


Lew 


S99 






47.26 


365M 


B 0457+24 


2.384 


S99 


8.8 


-0.77 


Lciv+V{CIY) 


Lew 


S99 






46.61 


365M 


B 2126-158 


3.270 


S99 


9.2 


-0.56 


iciv+V(CIV) 


Lew 


S99 






46.30 


333M 


B 2134+004 


1.936 


S99 


9.1 


-0.77 


iciv+V(CIV) 


Lew 


S99 






46.88 


365M 


OQ 172 


3.544 


098 


9.6 


-0.10 


Lh/3+V(H/3) 


Lhi3 


H03 


2.97 


H03 


47.61 


408M 


4C 12.50 


0.124* 


098 


7.8 




A/bh-o-* 




D06 


2.77 


G94 


44.62 


178M 



References: S89: Spencer et al. (1989); H08: Holt, J. et al. (2008); M97: Morganti et al. (1997); J05: Johnston et al. (2005); E04: Edwards & Tingay 
(2004); G06: Gallo et al. (2006); K06; Komossa et al. (2006); S03; Snellen et al. (2003); S99; Snellen et al. (1999); 098: O'Dea (1998); M06; McLure 
et al. (2006); J91; Jackson & Browne (1991); C97: Corbin (1997); S07: Sikora et al. (2007); W02; Woo & UiTy (2002); L96: Lawrence et al. (1996); 
H03; Hirst et al. (2003); C91: Corbin (1991); D03: D'Elia et al. (2003); G94: Gelderman & Whittle (1994); X08: Xie et al. (2008); M04: Marchesini et 
al. (2004); B90: Barthel, Tytler, & Thomson (1990); P06: Peng et al. (2006); H06: Holt et al. (2006); D06: Dasyra et al. (2006). 
Note: 

(a) The redshifts ai'e selected from NED^, and the redshifts with stars ai'e updated in Holt et al. 2008. 

(b) Methods for calculating the BH masses in column 4 as described in section 3.1. 

(c) Data used to calculate the bolometric luminosities in column 5 as described in section 3.2. 

(d) References for related data used in columns 6 and 7. 

(e) The BH mass is calculated from the relation between ionizing luminosity and BH mass (refer D'Elia et al. 2003 for more details). 
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Figure 1. The histogram distribution of redshift. 

for the classification as CSS/GPS sources. Columns (4)-(8) list the 
BH masses, Eddington ratios, methods for deriving the BH masses, 
methods for deriving bolometric luminosities, and references for 
related information in columns 6 and/or 7. Columns (9)-(10) give 
crjo III] and their references. Columns (11)-(12) list the jet power 
obtained from NEE0, and the waveband in which the power was 
calculated. 



3 METHODS 

3.1 BH mass estimation 

For the RM method, it takes a long time to simultaneously monitor 
the variability of the broad emission line and the continuum, and 

^ NASA/IPAC Extragalactic Database; |http://nedwww.ipac.caltech.edu | 



then obtain the BLR size. Fortunately, Kaspi et al. (2000) found an 
empirical correlation between the BLR size and the monochromatic 
luminosity at 5100A. Recently, Kaspi et al. (2005) reinvestigated 
the relationship between the characteristic BLR size, Rblh, and 
the H/3 emission line, X-ray, UV and optical continuum luminosi- 
ties. Several other empirical relations between i?BLR and emission 
line luminosities (e.g., Hq, C IV and IVlg II) have also been in- 
vestigated (e.g., Vestergaard 2002; Greene & Ho 2005a; Kong et 
al. 2006; IVlcGill et al. 2008). The BH mass then can be estimated 
with the assumption of virialized gas around the BH, 



Mbh = i?BLRl^'G~\ 



(1) 



where the quantity V can be estimated from the FWHM of the 
broad emission lines {V = /Vfwhm), and / = y/3/2 assuming an 
isotropic distribution of the BLR clouds (e.g., Wandel et al. 1999; 
Kaspi et al. 2000). 

The optical continuum in the radio loud (RL) AGN may be 
contaminated by the nonthermal emission from a jet: their optical 
luminosity and hence the viral BH mass may be over-estimated. 
In general, steep-spectrum radio sources tend to be oriented with 
jets pointing away from our line of sight, and the jet contribution at 
optical waveband may be negligible due to the weakness of beam- 
ing effects. Fanti et al. (1990) proposed that CSS sources were 
not strongly affected by Doppler boosting. Woo & Urry (2002) 
found that the BH masses derived from the optical luminosities 
and FWHM of the broad emission lines in RL AGN are consistent 
with those derived from other independent methods. Therefore, BH 
masses derived from broad line width and optical continuum/line 
luminosities for CSS sources should be credible. However, Baker & 
Hunstead (1995) found that the composite GPS quasar spectrum is 
clearly different from that of CSS quasars, but similar to that of flat- 
spectrum radio quasars (FSRQs). If this is true, the optical emission 
in GPS sources may be strongly contaminated by the synchrotron 
emission from the jets, as in FSRQs (e.g., Liu et al. 2006). There- 
fore, we use broad line luminosities to calculated the BH masses for 
GPS sources in our sample. We list the equations used to estimate 
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Figure 2. The correlations between the bolometric luminosities and broad line luminosities, Hf3 (left), Mg II (middle), and C IV (light) respectively for a 
sample of radio quiet AGN from Kaspi et al. 2005, where the bolometric luminosity Lbol = 9AL;^(5100A) is assumed (e.g., Kaspi et al. 2000). The solid 
lines are their best fittings. The optical continuum 5I00A luminosities and H/3 luminosities are selected from Kaspi et al. 2005, while Mg II, C IV luminosities 
are selected from Kong et al. 2006. 



BH mass from the broad emission line (Ha and H/3 from Greene 
& Ho 2005a; Mg II and C IV from Kong et al. 2006), 
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BH masses for some CSS sources, estimated from the optical lu- 
minosity and broad line width, are also selected from the literature 
directly if there is no line luminosity, but we correct the optical 
luminosity to our cosmology (see Table 1). 

The BH mass can also be derived from the A/bh — o"* relation 
(e.g., Tremaine et al. 2002) for nearby AGN (e.g., Gebhardt et al. 
2000b), 



Mf7 



0' 



(6) 



with an intrinsic scatter ~ 0.3 dex . 

Alternatively, the BH mass can be estimated from the galactic 
host bulge luminosity. McLure & Dunlop (2002) found that the BH 
mass is tightly correlated with the luminosity of the host galaxy 
for a sample of both active and inactive galaxies, with a scatter 
~ 0.4 dex. We use the relation between the host galaxy absolute 
magnitude Mr at Ji-band and BH mass proposed by McLure et al. 
(2004), 



logio(Mbh/M0) = -0.5Mh - 2.74, 



(7) 



to estimate the central BH masses for the sources without data for 
broad emission lines or bulge stellar velocity dispersion. 



3.2 Bolometric luminosity estimation 

The Bolometric luminosity of AGN is sometimes approximated 
from the optical luminosity since integration of the spectral en- 
ergy distribution (SED) is usually hampered by the lack of multi- 
wavelength observational data that spans many decades in wave- 
length. It is found that the bolometric luminosity is roughly 10 
times the optical luminosity, Lboi — lOLopt, where bolometric 
luminosity was integrated for the SEDs of AGN (e.g., Kaspi et al. 
2000; Woo & Urry 2002; see also Vasudevan & Fabian 2007). In 
this work, we adopt Lboi ~ OLsioo, where Lsioo = ALa(5100A) 
(Kaspi et al. 2000). However, the monochromatic optical luminos- 
ity cannot be used to derive the bolometric luminosity if beaming is 
important in radio loud AGN (e.g., BL Lac, FSRQs and maybe also 
GPS sources in our work). We know that the broad line emission 
is illuminated by the ionizing luminosity from the central engine, 
and, therefore, it can be used as a good indicator of thermal opti- 
cal emission. We can estimate the thermal optical emission using 
the broad emission lines for the sources in which the optical emis- 
sion may be contaminated by the jet (BL Lacs, FSRQs, and/or GPS 
sources). To do so, we fit the bolometric luminosities (assuming 
Lhoi = 9L5100) and the broad line luminosities (H/3, Mg II and C 
rV) for a sample of radio quiet AGN in Kaspi et al. (2005), where 
the optical emission is believed to be free from contamination by 
the nonthermal synchrotron emission (the jet is weak, if present). 
The H/3 luminosities and 5100A monochromatic luminosites are 
selected from Kaspi et al. (2005) for 35 AGN, of which 27 sources 
with Mg II luminosities, and 32 sources with C IV luminosities are 
selected from Kong et al. (2006). The relation between the bolo- 
metric luminosities and different broad emission lines (H/3, Mg II 
and C IV) are shown in Figure 2, and their linear regressions are 

log Lboi = 0.94 ± 0.03 log Lhp( + 5.39 ±1.14, (8) 

log Lboi = 0.74 ± 0.07 log LMgii + 13.06 ± 3.02, (9) 

log Lboi = 0.74 ± 0.06 log Lciv + 13.31 ± 2.61. (10) 

With Eqs. 8-10, we can calculate the bolometric luminosities for 
the radio loud AGN from their broad emission line luminosities, 
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especially for sources that are probably significantly contaminated 
by the beamed emission from the relativistic jets. We note that 
our bolometric luminosities derived from the H/3 luminosities are 
roughly 3 times larger than that derived using the popular rela- 
tion between the bolometric luminosity and the BLR luminosity 
Lboi — IOLblr, where Lblr ^ 25.26Lh/3 (e.g., Celotti et al. 
1997; Liu et al. 2006 and references therein). 



3.3 Jet power estimation 

Jet power is a fundamental parameter reflecting the energy trans- 
port to large spatial scales from the central engine by the radio jet. 
The jet power can be estimated from the low frequency radio emis- 
sion, which is usually emitted from the extended optically thin ra- 
dio lobes, and is free from Doppler beaming due to the low bulk 
velocity of the lobe motion (e.g., Willott et al. 1999). In this work, 
we use the formula proposed by Willott et al. (1999), 

^/^L15l''/^(W), 



Q. 



3 X lO^V 



(11) 



where L151 is the radio luminosity from the lobes measured at 151 
MHz in units of 10^**W Hz"^sr"^ Willott et al. (1999) have ar- 
gued that the normalization is uncertain and introduced the factor 
/ (1 ^ / ^ 20) to account for these uncertainties. Blundell & 
Rawlings (2000) argued that / ~ 10 is a likely consequence of 
the evolution of magnetic field strengths as radio sources evolve. In 
this work, we adopt / = 10, and search the low frequency radio 
data from NED for CSS/GPS sources. The extended emission was 
A'— corrected to 151 MHz assuming a — 1 {F^ oc v~"), and the 
jet power was then calculated from Eq. 1 1 (see Table 1). 



4 RESULTS 

4.1 BH mass, Eddington ratio and jet power 

The BH masses of 65 compact radio galaxies (38 CSS-l-27 GPS 
objects) are calculated in this work, or collected from the literature, 
using several different methods as described in Sect. 3.1. We find 
that the BH masses vary by several orders from 10'''^-10^°'^Mq, 
for both CSS and GPS sources. The histogram distributions of the 
BH masses are shown in Figure 3, the dotted-line, dashed-line and 
solid-line are for CSS sources, GPS sources, and the whole sample 
respectively. We find that the average logarithmic BH masses are 
<log i\fbh>=8.31, 8.23, 8.28 for 38 CSS sources, 27 GPS sources, 
and the whole sample respectively. 

The correlation between the bolometric luminosities and 
broad emission line luminosities (H/?, Mg II and C IV) for a sample 
of radio quiet AGN are explored in this work (see Figure 2 and Sect. 

3.2 for more details). We then calculate the bolometric luminosi- 
ties for all GPS sources, and some of CSS sources from their broad 
emission lines with the above empirical correlations (Eqs. 8-10), 
except for those with bolometric luminosities derived from their 
SED fittings. The Eddington ratios are calculated with I/boi/-i'Edd, 
where LEdd = 1-38 x lO^^Mbh/MQ. We find that young radio 
galaxies (GPS/CSS sources) have relatively high Eddington ratios, 
10~'^< iboi/iEdd< several (see Figure 4). The average logarith- 
mic Eddington ratios are <log Lboi/iEdd>=-0.55, -0.58 and -0.56 
for the 32 CSS sources, the 19 GPS sources, and their sum, respec- 
tively. 

It is found that the jet power ranges from ~10*'^ to 
~10*^erg s"^ for both the CSS and GPS sources. The histogram 
distributions of the jet power for the CSS/GPS sources are shown in 
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Figure 3. The histogram distribution of the black hole masses. 
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Figure 4. The histogram distribution of the Eddington ratios. 



Figure 5. The average logarithmic jet powers are <log Qjct>=46.1 
and 45.3 for the 37 CSS and the 26 GPS sources respectively. The 
average jet power in the CSS sources is larger than that of the GPS 
sources based on our sample. 



4.2 [O III] kinematics and accretion/jet properties 

Figure 6 plots BH mass, Mbh, versus velocity dispersion, crjo m], 
derived from FWHM(0 III), which includes 23 CSS (squares) and 
5 GPS sources (triangles). The solid line is the relation of Tremaine 
et al. (2002) obtained for 31 nearby inactive galaxies. It is clear 
that the CSS sources systematically deviate from the Tremaine et 
al. relation. Two of the five GPS sources also deviate substantially 
(OQ 172 and 4C 12.50). 

We calculate the deviation Act = log oyo ni] — log o"[prcd] 
for our sample with the measured ayo nij, where (J[pred\ is calcu- 
lated from the Tremaine et al. relation (Eq. 6) with the estimated 
BH mass. Figure 7a and 7b plot Act versus jet power, Qjct, and 
bolometric luminosity, iboi, respectively. Using the Spearman rank 
correlation analysis, there are no significant correlations between 
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Figure 5. The histogram distribution of the jet power. 
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Figure 6. The velocity dispersion ayQ ui] vs. BH mass Af^ij. The solid 
line is the — ct* relation of Tremaine et al. 2002 for 3 1 nearby inactive 
galaxies, and the dotted lines are their intrinsic scatter ±0.3 dex. 



A(7 and Qjct {r — 0.12 and Pnuii=0.51) or between Act and Lboi 
(r = 0.11 and Pnun=0.59). 

Figure 8 plots Act versus Eddington ratio I/boi/iEdd for the 
compact radio galaxies (20 CSS + 3 GPS sources). We find a signif- 
icant correlation between Aa and Lboi/iEdd for these 23 young 
radio sources, with the Spearman correlation coefficient r — 0.53 
and probability Pnuii — 1 x 10""^. For comparison with the ra- 
dio quiet AGN (the jet is weak, if present), we also calculate 
Act and the Eddington ratio I/boi/^Edd for 31 NLSls (circles) 
and 49 QSOs (crosses). The NLSls are selected from Grupe & 
Mathur (2004) and QSOs are selected from Shields et al. (2003), 
where BH masses are calculated from L5100 and the broad line 
width, and bolometric luminosities are calculated from Lsioo- We 
find that the relation, AfT-Lboi/iEdd, in the young radio galaxies 
(GPS/CSS) is similar to that of radio quiet AGN, and the sources 
with iboi/iEdd~l have the largest deviations (see Figure 8). 
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Figure 7. Jet power Qjct vs. Act (top), and accretion power Lbol vs. Act 
(bottom), where Act = logCT[o iii] - logCTfp^od] ■ 



5 DISCUSSION 

In this paper we investigate the central engine properties for a sam- 
ple of young radio galaxies (CSS/GPS sources), i.e., BH masses, 
Eddington ratios, jet powers, and then explore the possible relation 
between kinematics of [O III] narrow line and accretion/jet activi- 
ties. 

The BH masses of the young radio galaxies are estimated 
from several empirical relations in this work, as described in Sect. 
3.1, and some BH masses are selected from the literature after 
correcting to our cosmology, if the luminosities are used in de- 
riving BH masses. We find that most of the BH masses are in 
the range of IO^-IO^^Mq for both CSS and GPS sources (Fig- 
ure 3), which are typical BH masses for normal AGN. The mean 
BH mass of the 38 CSS sources is similar to that of the 27 GPS 
sources in our sample, with no evident difference between them. 
More complete sample is desired to further test this issue. However, 
the mean BH mass of young radio galaxies (CSS/GPS sources, 
< log Mbh >— 8.3) is systematically less than that of radio loud 
QSOs « log A/bh >~9. 1, e.g., Marchesini et al. 2004) or low red- 
shift radio galaxies « log A/bh >^S.9, e.g., Bettoni et al. 2003). 
One of the possible physical reasons is that the BHs in young radio 
galaxies are still growing rapidly. 

The bolometric luminosities for young radio galaxies are es- 
timated from SED fitting, monochromatic optical luminosities or 
broad line luminosities (see Sect. 3.2 for details). The bolometric 
luminosities of some CSS/GPS sources are not estimated due to a 



© RAS, MNRAS 000,[T]-?? 



8 Q.W.Wu 



0.5 



b 



-0.5 



n \ \ I \ \ 
■ CSS 

GPS 
X RQ QSO 
O RQ NLSl 



_L 



O 



_J \ I 1 \ \ \ I I \ L_ 



-2-10 1 
log U 



Figure 8. The relation between the Eddington ratio ibol/^Edd and Act. 

lack of information on their classification (e.g., type I or II) where 
optical luminosities and then bolometric luminosities may be un- 
derestimated if they are obscured by the putative torus. Normally, 
steep-spectrum radio sources do not suffer from strong contami- 
nation from the Doppler-boosted emission from the jet, due to the 
possible large inclination angles. However, there is some evidence 
from radio variability that at least some of the GPS sources are 
Doppler-boosted (e.g., O'Dea 1998, and references therein). The 
optical composite spectrum of some GPS sources is also clearly 
different from that of CSS sources, which is similar to that of FS- 
RQs and may be contaminated by the nonthermal emission from 
jets (e.g.. Baker & Hunstead 1995). Therefore, the bolometric lu- 
minosities of all GPS sources are calculated from the broad line lu- 
minosities (if present) using our empirical correlations constrained 
from a sample of radio quiet AGN (Eqs. 8-10). We estimate the 
Eddington ratios, Lboi/^'Edd, from estimated BH masses and bolo- 
metric luminosities, and find that all CSS/GPS sources are high Ed- 
dington ratio systems (10~'^<Z/boi/iEdd<several, see Figure 4). 
The mean Eddington ratios are < log Lboi/i'Edd >=-0.55 and - 
0.58 for CSS and GPS sources respectively, which are higher than 
that of broad line radio galaxies < logLboi/iEdd >=-1.9I or ra- 
dio loud QSOs (< log Lboi/i'Edd >=-l.ll, e.g., Marchesini et al. 
2004), but similar to those of NLSls (< log Lboi/I/Edd >~-0.8, 
e.g., Bian et al. 2008). Our results suggest that these young radio 
galaxies are also in the early stage of the accretion activity and that 
the BHs are still growing rapidly. Czemy et al. (2009) proposed that 
the short-lived young radio galaxy may be due to the intermittent 
activity caused by the radiation pressure instability in the accretion 
disk, where the instability occurred when the accretion rate (in the 
Eddington unit) exceeds a few percent. It should be noted that most 
of the young radio galaxies in our sample are high Eddington ra- 
tio sources (see Figure 4), which may have suffered the radiation 
pressure instabilities. If this is true, the radiation pressure instabil- 
ity will offer an interesting physical mechanism for the possibly 
short-lived young radio galaxies (Wu 2009, in preparation). 

Kawakatu et al. (2009) found that young radio loud AGN fa- 
vor an accretion disc without the big blue bump of the standard 
accretion disc, based on photoionization modeling of various nar- 
row emission lines, and suggested that the accretion mode in these 
young radio galaxies may be though advection dominated accre- 
tion flow (ADAF, see Narayan & McClintock for a recent review 
and references therein). It is widely beheved that optically thin 



ADAF exists only for Eddington ratios less than a critical value 
(Lboi/iEdd)c — 0-01 (e.g., Narayan & Yi 1995), and the accre- 
tion mode should transit to an optically thick standard accretion 
disc when the Eddington ratio larger than this critical value (e.g., 
Wu & Gu 2008 and references therein). However, our results show 
that nearly all the young radio galaxies in our sample have Edding- 
ton ratios larger than this critical value, and, therefore, it does not 
support the ADAF scenario. The difference in the narrow line ratios 
between young radio loud AGN and radio quiet Seyfert 2 galaxies 
in Kawakatu et al. (2009) may be caused by other factors. 

The width of the [O III] emission line has been used as a proxy 
for stellar velocity dispersion to derive the BH masses in AGN, 
even through the scatter is large (e.g.. Nelson 2000; Boroson 2003; 
Shields et al. 2003; Bonning et al. 2005). However, previous stud- 
ies of NLSl galaxies, employing different samples and methods, 
led to partially conflicting results in that the NLS 1 galaxies are sys- 
tematically off the Mhh — o"* relation if cr, = ayo nij is assumed 
(e.g., Bian & Zhao 2004; Grupe & Mathur 2004; Komosssa & Xu 

2007) . The [O III] width seems to be broader than gravitational ve- 
locities in these NLSls. It has been shown that CSS sources also 
tend to have rather broad [O III] emission lines (e.g., Gelderman 
& Whittle 1994; Nelson & Whittle 1996; Tadhunter et al. 2001; 
O'Dea et al. 2002; Holt et al. 2008). However, a direct comparison 
of [O III] width with the central BH mass has not yet been made 
so far. We estimate the BH masses for a sample of young radio 
galaxies (CSS/GPS sources) in this work, and find that most of the 
sources lie systematically ojf the Afbh — o"* relation, when using 
(T» = ayo HI], as found for NLSls (see Figure 6). A difference 
is that more than 50% of BH masses in young radio galaxies are 
larger than IO^Mq and many of them have broad line widths of 
> 2000 km s~^, while most BH masses of NLSls are less than 
IO^Mq with broad line widths of < 2000 km s"^ (e.g., Grupe 
& Mathur 2004). Therefore, the BH mass is not the key parame- 
ter causing the deviation of [O III] emission line width. We note 
that the extreme [O III] kinematic components, as strong wings 
with blueshifts up to ~2000 km/s, are found in both radio quiet 
NLSls (e.g., Aoki et al. 2005; Bian et al. 2005; Komossa et al. 

2008) and young radio galaxies (e.g., Tadhunter et al. 2001; Holt et 
al. 2008), and these strong blue wings are often not well included in 
the FWHM fittings. However, the sources with extreme kinematics 
are still rare (e.g., Bian et al. 2005), which will not change our main 
conclusion. 

In the standard model, the narrow emission lines are excited 
through photoionization by the UV continuum produced by the 
accretion disc. The situation is more complicated in radio galax- 
ies, since shocks caused by the interaction of jet and interstel- 
lar medium may also be important for the morphology, kine- 
matics, and excitation of the narrow emission lines. Inskip et al. 
(2002) investigated eight 6C radio galaxies at redshift 2 ~ 1 and 
concluded that combination of accretion disc photoionization and 
shock photoionization provide the best explanation of their spec- 
tra. Furthermore, by modeling the ratios of different narrow emis- 
sion lines with shock and AGN photoionization models, Moy & 
Rocca-Volmerange (2002) suggested that the shock ionization may 
be dominant for radio sources with intermediate radio sizes D (2 
kpc<Z)<150 kpc, e.g., CSS sources), while AGN photoionization 
is dominant by the most compact radio sources (D<2 kpc, e.g., 
GPS sources) and the extended sources with Z)>150 kpc. We note 
that there are two GPS sources (OQ 172 and 4C 12.50) with ra- 
dio sizes of less one kpc, which also show large deviations in their 
[O III] width (Figure 6). More compact GPS sources are needed to 
further test this issue. 
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The jet-cloud interaction is expected to produce extreme kine- 
matics of narrow emission lines, and the jet-driven outflows pro- 
vide a convenient mechanism to explain the highly broadened line 
profile, large velocity shift, and strong alignment between radio jet 
and optical line in young radio galaxies (e.g., de Vries et al. 1997, 
1999; Axon et al. 2000; Tadhunter et al. 2001; Holt et al. 2003, 
2006, 2008; Labiano et al. 2005). We investigate the possible re- 
lation between the deviation of [O III] line width Aa and the jet 
power Qjot of young radio galaxies. However, there is no signifi- 
cant correlation between them based on our sample. Therefore, the 
jet power should be not the only parameter caused the extreme [O 
III] kinematics if the jet-driven scenario is indeed important. 

The Eddington ratio, which is proportional to the accretion 
rate per unit BH mass, characterizes the extent to which radiation 
pressure competes with gravity in the nucleus, and may also play 
an important role in shaping the kinematics of the [O III] emission 
line apart from the primary driver of gravitational potential of the 
bulge (e.g., Greene & Ho 2005b). After removing the blue wings, 
the core component of oyo nij is still a good proxy of stellar ve- 
locity dispersion a* (e.g., Greene & Ho 2005b; Komossa & Xu 
2007). Blue wings are always found in NLSls with high Eddington 
ratios (~ 1, e.g., Komossa et al. 2008), which may result from the 
winds/outflows driven by the strong radiation pressure in these high 
Eddington ratio sources. It is interesting to note that most of young 
radio galaxies also have high Eddington ratios, which are similar 
to those of NLSls (Figure 3). We find that the deviations of the [O 
III] emission line in these young radio galaxies, Atj, are also well 
correlated to the Eddington ratios (see Figure 8). In particular, the 
young radio sources with Eddington ratios Lboi/iEdd ~ 1 have 
the largest deviations, which are also similar to those of radio quiet 
NLSl galaxies (see Figure 8). The similarity in [O III] width devi- 
ations between radio loud CSS/GPS sources and radio quiet AGN 
(QSOs and NLSls) suggests that the accretion activities may play 
a more important role than the jet activities in shaping the kinemat- 
ics of [O III] narrow emission line, since the jet is absent or weak 
in radio quiet AGN. The alignment of the narrow emission-line gas 
with the radio jet is always used to as an evidence for the shock 
scenario (e.g., de Vries et al. 1997; de Vries et al. 1999; Axon et 
al. 2000; Privon et al. 2008). However, it should be also natural for 
the alignment of radio jet and [O III] emission line if the outflow is 
driven by the radiation pressure caused by accretion activities. 



6 SUMMARY 

Using a sample of 65 young radio galaxies (38 CSS -I- 27 GPS 
sources) selected from the literature, we estimate BH masses, Ed- 
dington ratios, and investigate the possible physical reasons for de- 
viations of width of [O III] narrow emission line. The main results 
can be summarized as follows: 

(1) The BH masses in young radio galaxies (CSS/GPS 
sources) range from IO^'^-IO^'^'^Mq, which is the typical value 
in AGN (see also Gu et al. 2009), and the mean value is 
<log A/bh>— 8.3 which is less than that of radio loud QSOs or 
low redshift radio galaxies (< log Afbh >— 9.0, e.g., Marchesini et 
al. 2004; Bettoni et al. 2003). 

(2) Most of the CSS/GPS sources have relatively high Edding- 
ton ratios, with mean value <log I/boi/i'Edd>=-0.56, which are 
similar to those of NLSls. This result supports the idea that the 
young radio galaxies are not only in the early stage of jet activities 
but also in an early stage of accretion activities. 

(3) The A/bh — t7[o mj relation in young radio galaxies sys- 
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tematically deviates that defined by nearby inactive galaxies, and 
reminiscent of NLSls. We find that the deviation of the [O III] 
width Act is not related to the jet power or accretion power, but 
is related to the Eddington ratio (Figures 6, 7, and 8). There is no 
evident difference between the radio loud CSS/GPS sources and ra- 
dio quiet QSOs/NLSls in the Acr-Lboi/iEdd relation, which may 
suggest that accretion activities may still play an important role in 
shaping the kinematics of the [O III] emission line in these young 
radio galaxies, since that radio jet is weak or absent in radio quiet 
AGN. 
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